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Abstract 

We study the Electric Dipole moment of b quark in the general two Higgs Doublet 
model (model III) and three Higgs Doublet model with 0{2) symmetry in the Higgs 
sector. We analyse the dependency of this quantity to the new phase coming from the 
complex Yukawa couplings and masses of charged and neutral Higgs bosons. We see that 
the Electric Dipole moment of b quark is at the order of 10 -20 e cm, which is an extremely 
large value compared to one calculated in the SM and also two Higgs Doublet model 
(model II) with real Yukawa couplings. 
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1 Introduction 



The study of CP violating effects provides comprehensive informations in the determination of 
free parameters of the various theoretical models. Non-zero Electric Dipole Moment (EDM) 
for elementary particles is the sign of such violation. Neutron EDM has a special interest 
and the experimental upper bound has been found as < 1.1 x 10~ 25 e cm ]]J. EDM of 
electron and muon have been measured experimentally as d e = (— 2.7 ± 8.3) x 10~ 27 e cm || 
and = (3.7 ± 3.4) x lCT 19 e cm ||. These measurements can also give powerful clues about 
the internal structure of the particles, if it exists. 

The source of CP violation in the SM is complex Cabbibo-Kobayashi- Maskawa (CKM) 
matrix elements. The quark EDM vanishes at one loop order in the SM, since moduli of matrix 
element is involved in the relevant expression. Further, it also vanishes at two loop order after 
sum over internal flavours @, [5|. When QCD corrections are taken into account, nonzero EDM 
exists §. However, EDM for quarks is very small in the SM and need to be enhanced with the 
insertion of new models beyond. In the literature, quark EDM is calculated in the multi Higgs 
doublet models [7], |8|, Q. In ||, EDM due to the neutral Higgs boson effects is obtained in the 
two Higgs doublet model (2HDM) and || discusses the necessity of more scalar fields than just 
two Higgs doublets for non-zero EDM when only the charged Higgs boson effects are taken into 



account. In JT0[, EDM of 6-quark in the case of Z-boson and photon are calculated in the 2HDM 
and models with three and more Higgs doublets. EDM of b-quark for Z-boson is predicted in 
the range 10~ 21 — lCT 20 e cm, following the scenario where CP violation may only come from 
the neutral Higgs sector. Further, 6-quark EDM is obtained in the range 1CT 23 — lCT 22 e — cm 



when the CP violating effect comes from the charged sector. |T| is devoted to the study of 
t-quark EDM, which is predicted at the order of the magnitude of 10 _20 e cm, arises from the 
neutral Higgs sector in the framework of the multi Higgs doublet model. In [TJJ], quark EDM is 
calculated in the 2HDM if the CP violating effects come from CKM matrix elements and it is 
thought that H ± particles also mediate CP violation besides W bosons, however, at the two 
loop order these new contributions vanish. The electric and weak electric dipole form factors 
for heavy fermions in a general two Higgs doublet model is studied in |13[ and it is concluded 
that the enhancement of three orders of magnitude in the electric dipole form factor of the b 
quark with respect to the prediction of 2HDM I and II is possible 

In our work, we study EDM of 6-quark in the general 2HDM (model III) and the general 
3HDM with 0(2) symmetry in the Higgs sector (3HDM(02)) IQ- 111 this case, CP violation 
comes from complex Yukawa couplings and even at one loop, it is possible to get extremely large 
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EDM, at the order of 10~ 20 e cm. Since the theoretical values of quark EDM's are negligible in 
the SM and also in the model II 2HDM with real Yukawa couplings, model III and 3HDM(0 2 ) 
results can give important information about the new parameters, namely masses of charged 
and neutral Higgs particles, Yukawa couplings, beyond the SM. Note that, non-zero EDM can 
be obtained due to the charged Higgs boson effects even in the 2HDM in this case. 

The paper is organized as follows: In Section 2, we present EDM for b quark in the framework 
of model III and 3HDM(0 2 ). Section 3 is devoted to discussion and our conclusions. 

2 Electric Dipole moment of b quark in the general two 
Higgs Doublet and three Higgs Doublet models 

In this section, we calculate 6-quark EDM in the model III and extend our results to the general 
3HDM. Since nonvanishing EDM is the sign for the CP violation, we assume that there exist 
complex Yukawa couplings which are possible sources of such violations. We start with the 
Yukawa interaction in the model III 

£y = VijQiL^iUjR + VijQiL^iDjR + ^Q iL faU jR + gjQiL<hDjR + h.c. , (1) 

where L and R denote chiral projections L(R) = 1/2(1 =f 75), fa for % — 1, 2, are the two scalar 
doublets, QiL are left handed quark doublets, Ujn(Dj R ) are right handed up (down) quark 
singlets, with family indices i,j. The Yukawa matrices r^- D and £^- D have in general complex 
entries. Here fa and fa are chosen as 



01 = 71 



\ , ( V2x 
v + H° J + { t X ° 



^2 ( Ih + iH 2 ) ' (2) 



with the vacuum expectation values, 



<01>= 7! (°) ;< ^ 2> = ' ( 3 ) 

to collect SM particles in the first doublet and new particles in the second one. The Flavor 
Changing (FC) part of the interaction can be written as 

£y,FC = &QiL<foU jR + QQiLfaDjR + h.c. , (4) 

where the couplings £ U ' D for the FC charged interactions are 

€ch = £n V C km , 

£ch = v ckm £n > ( 5 ) 



and is denned by the expression (for more details see |15| ) 

Note that the index "N" in £jy D denotes the word "neutral". 
The effective EDM interaction for 6-quark is defined as 

C EDM = tdb l5 a^bF^ , (7) 

where is the electromagnetic field tensor and "d" is EDM of 5-quark. Here, "d" is a real 
number by hermiticity. In model III, the charged H ± and neutral Higgs bosons can 
induce CP violating interactions which can create EDM at loop level. Note that, we take Hi 
and H 2 as the mass eigenstates h° and A respectively and do not take the mixing of two 
CP-even neutral bosons Ho and ho into account, since no mixing occurs at tree level for our 
choice of Higgs doublets (eq. (0)). We present the necessary 1-loop diagrams due to charged 
and neutral Higgs particles in Figs. [I] and |[ Since, in the on-shell renormalization scheme, the 
self energy J2(p) can De written as 

J2(p) = (p - m b )J2(p)(P - m h ) , (8) 

diagrams a, b in Figs. [I] and |2| vanish when 6-quark is on-shell. However the vertex diagrams c, 
d in Fig. [I] and c in Fig. |2| give non-zero contributions. The most general vertex operator for 
on-shell 6-quark can be written as 

r M = F 1 (g 2 ) 7M + F 2 (g 2 )a^g^ 

+ F 3 (q 2 )a^ 5 q» (9) 

where q v is photon 4- vector and q 2 dependent form factors Fi(q 2 ) and F2(q 2 ) are proportional 
to the charge and anamolous magnetic moment of 6-quark respectively. CP violated interaction 
exists with the non-zero value of F^(q 2 ) and it is proportional to EDM of 6-quark. By extracting 
the CP violating part of the vertex, 6-quark EDM "d" (see eq.([?|)) is calculated as a sum of 
contributions coming from charged and neutral Higgs bosons, 

d = d H± +d h ° + d A ° , (10) 

where d H± , d h ° and d A ° are 

H ± 4 Gf e 1 - v - D 2 y((-l + Q t (-3 + y) - y)(y - 1) + 2 (Q t + y)lny) 

d = 7^ oTTT — ^N,tt lm {$N,bb) \ v tb\ 



V2 32tt z m t ' (y - 1 
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V2 167r 2 m b J n ( ri _ 4) 2 2 

d A ° = -d h \ ri ^r 2 ), (11) 



-. Here Q b and Qt are charges of b and t quarks 



with r\ = mlo/ml, r 2 = m 2 A0 /m 2 and y = ■ 
respectively and t^ff = J • in ec i- (0) we take into account only internal t-quark 
contribution for charged Higgs and internal 6-quark contribution for neutral Higgs interactions 



since we assume that the Yukawa couplings i = u,c, Cjv,bj ? J : = d, s and £jy tc are negligible 
compared to ^ tt and ^Nbb ( see ED- Further, we choose ^tt Tes ^ an d £/v&6 complex, 



ED 
SAT, 66 



l£.W,6bl € 



(12) 



Finally, using the parametrization eq.(|T^), we get the EDM of 6-quark in model III as 



d = ^ , \^ bb \sin9{— & tt \V tb \ 2 



V2 32tt 2 

+ 2—— l^N bbl COS I 

r 2 (r 2 - 2) 



y ((-1 + Q t (-3 + y) - y)(y - 1) + 2 (Q t + y) Iny) 

(y-i) 8 



J_ =7 

ri (ri - 2) 0=1 - V r i - 4 , 1 . 

/ m- o o ri /nri 

yVi (n - 4) ^ ^ 



yV 2 (r 2 - 4) 



, v 77 ^ - V^2 - 4 1 
m— -r 2 tnrz)} ■ 



(13) 



Now, we would like to extend our result to the general 3if-DM(0 2 ) [[14]]. The new general 
Yukawa interaction is 



C 



(14) 



where 0, for z = 1,2, 3, are three scalar doublets and 7}^°, , p^- D are the Yukawa matrices 
having complex entries, in general. With the choice 



1 

71 



i 



and the vacuum expectation values, 





v + H° 



( V2r 



(15) 



J_ / V2F+ 



< 01 >: 



- 

V2 V v 



< 



>= ; < </>, >= 



(16) 



the information about new physics beyond the SM is carried by the last two doublets 2 and 
03, while the first doublet <j>\ describes only the SM part . The Yukawa interaction for the 
Flavor Changing (FC) part is 



£y,fc = QiLfoUjR + ^QiL&DjR + PijQiLfoUjR + PijQiLfoDjR + h.c. 
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(17) 



where, the charged couplings ^ D and p^ D are 



CKM 



P^/i = Pn Vckm , 

Pch = V C km Pn , (18) 

and 

<:U{D) _ (U U(D),-1MD) V U(D) 
SJV — \ V R(L) ) S L(R) > 

PT ] = (V^)-V^V^ . (19) 

In the 3HDM, there are additional charged Higgs particles ,-F ± , and neutral Higgs bosons h'°, 
A' (see ||14|| ). These particles bring new part to EDM of 6-quark and by taking only couplings 
sNtt ^Nbbi Pnu an d PNbb account, it can be written as 

d = d H± + d h ° + d A ° + d F± + d h '° + d A '° , (20) 
where d F± , d h '° and d A '° are 

,f± 4 Gp e 1 _ v D , T/ 2 y ((-1 + Q t (-3 + y') - y%/ - 1) + 2 (Q f + y') jwyQ 

" ^ 32vr^ P ^ Jm(p ^ )l ^ 1 {^Tf 



V2 16vH m 6 / r / ( r / _ 4) 2 2 

d* = -dVi-ri), (21) 

with r' x = m\,o/m b , r 2 = m 2 A , /ml and y' = The expressions for d H± , g^ and c?" 40 are 

defined in eq.(|ll|). Further, by introducing 0(2) symmetry in the Higgs sector, the masses of 
new Higgs bosons F , h'° and A' become the same as the masses of model III Higgs bosons, 
H , h° and A respectively (see fl4}l). Therefore, the final result for EDM of 6-quark in the 
3HDM(0 2 ) is 

d = ^ ^ {^(S^M^ +PN, tt Im(PN, bb )) \V»\~ 



y((-l + Q f (-3 + y)-y)(y - 1) + 2 (Q f + y)lny) 
+ 2 ^ (Im(^ bb ) Re(i° bb ) + Im(p N , bb ) Re(p N>bb )) 



ri(ri-2) v/rT - aAi - 4 , 1 , 
( / = In- o + o n Zrari 

yjn (n - 4) 2 2 

7-2(7-2-2) ^ - y/r 2 - 4 1 

y . . ^ ^r 2 lnr 2 )}. (22) 

^r 2 (r 2 - 4) 2 2 
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New 0(2) symmetry also permits us to parametrize the Yukawa matrices £ U ( D ^ and jFW as 

El 

^U(D) =e U(D) cose ^ 

p u = e u sin 9 , 

p D = ie D sin9, (23) 
where are real matrices satisfy the equation 

Here T denotes transpose operation. In eq. (|2"B|), we take p D complex to carry all CP violating 
effects on the third Higgs scalar. Using the parametrization in eq. (|23|) , it is easy to see that 
only charged Higgs bosons contribute to EDM of 6-quark but not the neutral Higgs ones, since 
f>Nbb 1S a pure imaginary number. Finally, 6-quark EDM reads as 

, 4GV e l |T/|2 _c/ _ D ■ 2 ^((-l + Qt(-3 + y)- y)(y -l) + 2(Q t + y)lny) ln ^ 

3 Discussion 

In this section, we study dependencies of 6-quark EDM on the masses of charged and neutral 
Higgs bosons and the CP violating parameter 9. In the analysis, we use the CLEO measurement 

E3 

Br(B -> X s7 ) = (3.15 ± 0.35 ± 0.32) 10" 4 . (26) 

to find a constraint region for our free parameters. The procedure is to restrict the Wilson coef- 
ficient Cy , which is the effective coefficient of the operator 7 = j^s a a^ u (m b R + m s L)b a J-'v L/ 
(see 



18| and references therein), in the region 7 < \C e 7 ff \ < 0.439. Here r and lower 
limits were calculated using the CLEO measurement and all possible uncertainities in the cal- 
culation of This restriction allows us to find a region for the parameters ^N tt , C^b&> ^ 
in the model III and e^ tt , e® bb , 9 in the general 3HDM(02). In our numerical calculations, we 
also respect the constraint for the angle 9 due to the experimental upper limit of neutron elec- 
tric dipole moment, d n < 10 _25 e-cm, which leads to -^—I , m(^ tt ^* I f bb ) < 1.0 for M H ± rs 200 



GeV Pf in the model III and -±— {e u N tt e% bb ) sin 2 9 < 1.0 in 3HDM(0 2 ). Further, we take 



m t m b 

^^^N,tt c N,bbJ " L 

only and {^nu an< ^ ^Nbb) nonzero in the model III (3if-DM(0 2 )) and neglect all other 
couplings. Note that h° is assumed as the lighest Higgs boson in our calculations. 



6 



In Fig. ^ we plot EDM n d n with respect to the ratio R neu tr — for sin 9 = 0.5, m A o = 
80 GeV, m H ± = 400 GeV, £§ bb = 40 m b and \r tb \ = \ < 1 in the model III. Here "d" lies in 

' $N,bb 

the region bounded by solid lines for > and by dashed lines for < 0. It is observed 
that 6-quark EDM is strongly sensitive to the ratio R neu tr and this dependence increases with 
decreasing masses of neutral Higgs bosons. If the ratio R neu tr becomes small the considerable 
enhancement of d will be obtained and therefore the mass difference of h° and A should not 
be large. In the case of degenerate masses of h° and A , the contribution of the neutral Higgs 
sector to 6-quark EDM vanishes. 

Fig. H (^) is devoted to sin 9 dependence of "oT for m h o = 70 GeV, m A o = 80 GeV (m h o = 
m^o = 80 Gel^), m H ± = 400 GeV, £® bb = 40 m h and |r t &| < 1. If the value of sin 9 decreases, 
n d" becomes small as expected and the restricted region becomes narrower, for both C 7 e// > 
and C 7 e// < 0. "tf is positive for C 7 e// > 0, however, for C 7 e// < 0, it can also take 
negative values when m^o reaches to m A o . This is an intersting result which can be used in the 
determination of the sign of Gj. 

In Fig. H we plot "rf" with respect to the charged Higgs mass m#± for sin 9 = 0.5, 
m h o = 70 GeV, m A o = 80 Ge^, ^ bb = 40m b and \r tb \ < 1. This figure shows that "d" is 
weakly sensitive to the charged Higgs mass for m#± > 400 GeV, especially in the case 
< 0. The restriction region becomes narrower with increasing m#±. 

In 3HDM(02), the possible choice of the parametrization for Yukawa couplings (eq. (|23|) ) 
causes to vanish the contribution of the neutral Higgs sector to the 6-quark EDM (see eq. (p5|)) 
and therefore only the charged Higgs sector contributes. 

Fig. is devoted to sin 9 dependence of "d" for m H ± = 400 GeV, ^ bb = 40 m b and |r t &| < 1 
in 3if-DM(0 2 )- The behaviour of EDM is similar to the model III case with degenerate masses 
m^o and m^o, however, "d" is more sensitive to sin 9, since it is proportional to sin 2 9 (see 
eq. (^)). Further, "d" is not sensitive to the charged Higgs mass 777#± for m H ± > 400 GeV, 
especially in the case Cj^ < 0, similar to the model III. (Fig. ^) 

Now we would like to summarize the main points of our results: 

• It is interesting that EDM is generated by the one loop diagrams. This is due to the 
freedom to choose the Yukawa couplings as complex numbers in the models under consid- 
eration. Further, even in the model with two Higgs doublets, there exist a non-vanishing 
contribution to EDM due to the charged sector. 

• Since 6-quark EDM "d" is strongly sensitive to the ratio R neu t r , in model III, the mass 
difference of h° and A should not be large. For the 3HDM under consideration, there is 
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no need to restrict the neutral Higgs masses because they do not contribute to "d" for 
the given parametrization of Yukawa couplings. 

• If "d" is positive, can have both signs. However, if it is negative, C7 must be 
negative for both models. This is an important observation which is useful in the deter- 
mination of the sign of . 

• "d" is not sensitive to the mass of charged Higgs boson for its large values in both models. 

• EDM of "6" quark is at the order of ~ l(r 20 e cm in both model III and 3HDM{0 2 ) 
and its magnitude is larger compared to the results (~ 1CT 23 — lCT 22 e cm) existing in the 
literature. The neutral Higgs boson effects are strong in the model III and there is an 
enhancement, even one order (d ~ (1CT 19 ), if the masses of neutral Higgs bosons, m h o 
and m^o, are far from degeneracy. 

Therefore, the experimental investigation of the 6-quark EDM gives powerful informations 
about the physics beyond the SM. 
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Figure 1: One loop diagrams contribute to EDM of fe-quark due to H ± in the 2HDM. Wavy 
lines represent the electromagnetic field and dashed lines the H ± field. 
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Figure 2: The same as Fig. [I], but for neutral Higgs bosons h° and A . 
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m h Jm Ao 



Figure 3: ft-quark EDM "d" as a function of the ratio R neu tr = — b2 -, for m H ± = 400 GeV, 
sin9 = 0.5, £jy b6 = 40 m b and \r t b\ = \-nr L \ < 1, in the model III. Here d is restricted in the 

^N,bb 

region bounded by solid lines for Cj^ > and by dashed lines for Cj^ < 0. 
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Figure 4: 5-quark EDM "cf' as a function of sin# for m# = 400 GeV^, m^o = TOGeV^, m^o = 
80 GeV, = 40 rat, and | rt& | < 1, in the model III. Here d is restricted in the region bounded 
by solid lines for > and by dashed lines for C7 < 




Figure 5: The same as Fig|| but for m ft o = m^o = 80 GeV. 
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Figure 6: 5-quark EDM "d" as a function of for sin 6 = 0.5, m h o = 70 GeV, m^o = 80 GeV, 
£.Nbb = 40 m(, and |r#| < 1, in the model III. Here d is restricted in the region bounded by solid 
lines for Cy > and by dashed lines for < 0. 




Figure 7: The same as Fig. but in 3if_DM(0 2 ). 
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Figure 8: The same as Fig. |, but in 3HDM(0 2 ). 
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